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ABSTRACT 

We use the current weak lensing data to constrain the equation of state of dark energy w and 
the total mass of massive neutrinos yjm„. The constraint on w would be weak if only the current 
weak lensing data are used. With the addition of other observational data such as the type la 
supernovae, baryon acoustic oscillation, and the high rcdshift Hubble parameter data H(z) derived 
from relative galaxy ages to break the degeneracy, the result is significantly improved. For the pure 
uiCDM model without massive neutrinos, we find w — — I.OOIq'^. For the wCDM model with the 
massive neutrino component, we show that the constraint on w is almost unchanged, there is very 
little degeneracy between w and ^m„. After marginalizing over other parameters, we obtain the 
probability distribution function of ^ m„, and find that the upper limit is m„ < 0.8 eV at 95.5% 
confidence level for the combined data sets. Our constraints of w and ^2 m v are both compatible and 
comparable with the constraints obtained from the WMAP 5-year data. 

Subject headings: cosmology: theory — cosmological parameters — gravitational lensing — neutrinos 



1. INTRODUCTION 

Great progress is achieved in the study of modern 
cosmology in recent years. Thanks to the precision 
observations of the "distance indicator" type la su- 
pernovae (SN la) and the anisotropy of the cosmic 
microwave background (CMB), an accelerating Uni- 
verse with a non-baryonic da r k matter component is 
established dRiess et all 119981: iPerlmutter et all 119991 : 
Ide Bernardis et alj 12000 : ISpergel et alj 120031) . " 



cosmo logical probes 



Other 

such as the large scale struc- 
galaxy clusters (|Allen et al.l 



tures (|Tegmark et al.H200 

2008) , H ubble rate deriv ed from relative galaxy ages 
(RGA) ([Simon et ail 120051 ). provide further supporting 
evidence for this scenario, and help determine the model 
parameters with even better precision. However, the na- 
ture of the so-called dark energy which drives the acceler- 
ation of the Universe is still far beyond our understand- 
ing. The simple cosmological constant cold dark matter 
model (ACDM) remains to be the most popular, with 
an equation of state (EOS) w = p/p = — 1 for the dark 
energy. However, such a simple cosmological constant 
suffers from the fine tuning and co incidence problems 
([Weinberg! [19891 IZlatev et alj 1 19991) . Dynamical dark 
energy models typically have evolving EOS w — w(z), 
a wid e variety have been proposed and tested in litera- 
ture (ICaldwell 120021: IZhao et al.ll2005l ; iFeng et alj[2006t 
IXia et al.l 12000 : IZhao et all 120071 ) ~To eventually solve 
the dark energy problem, it is very important to develop 
additional techniques and use more observational data to 
test the dark energy models and measure other cosmo- 
logical parameters. 

Another fundamental problem in modern physics 
which we would like to draw particular attention to in 
this paper is the mass of neutrinos. The neutrino os- 
cillation experiments suggest that the mass differences 



between various types of neutrinos are Amf 2 ~ 8 x 10 -5 
eV 2 (from solar neutrino experiment) and Am 2 2 ~ 
2.2 x 10~ 3 eV 2 (from atmospheric neutrino experiment) 
(jMaltoni et all 120041 ). These results suggest that the 
masses of the neutrinos form a hierarchy of mi ~ 0, 
m 2 ~ Am solar and m 3 ~ Am at mosphcric, or alternative 
an inverted order. However, it is also possible for the 
three types of neutrinos have almost degenerate masses 
ni\ ~ fn2 ~ W3 3> Am a t m0 spheric- There are not abso- 
lute bounds on the neutrinos masses from oscillation ex- 
periments, while the cosmological observations can pro- 
vide us effective ways to infer the abundance and total 
mass of neutrinos, since the kinematics of n eutrinos af- 
fects the growth of cosmological structures (|Hannestadl 
2006). For the ACDM cosmology, upper bounds on the 
sum of the neutrino masses have reached the level of 
m u 5, 0-3 — 1 eV at 2a depending on the data set 
used (jHannestad et al.1l2006l and references therein). If 
the dark energy EOS is taken into account however, the 
upper boun d is shown to be relaxed to about 1.5 eV 
at 2a level (jHannestad 2005). This may indicate a de- 
generacy betweejithfi_j^utnno_jnass_a^^ energ y 
EOS w (|Hannestadl 120051 : iLesgourgues fe PastoH[2006f) . 
The inclusions of ba ryon acoustic oscillation (BAO) data 
( Goobar et al.l [20061) an d the weak gravitational lensing 
( Hannestad et al .1120061 ) may help break such degeneracy, 
as we shall proceed to show in the remaining part of this 
paper. 

The gravitational weak lensing (WL) effect generates 
small distortions (of the order 1%) on the images of 
distant galaxies as the light passed through inhomoge- 
nous matter distribution along the line of sight. This 
effect can be measured statistically to yield information 
on the density field and the geometry of the Universe. 
A great advantage of WL is that it relies on the total 
matter content of the Universe directly, thus the prob- 
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lem of galaxy-to-matter bias in the large structure sur- 
vey is avoided. The WL surveys provide a powerful and 
precise probe of the dynamical property of the Universe 
at redshift z < 3, and it is an independ ent technique 
to m easure the properties of dark energy ([Munshi et al.l 
2()()8 j. "Cosmic shear" is recently reported in galaxy 
WL surveys, which provides imp o rtant cosmo l ogica l 
implications (IWittman et alj [2000: iBacon et alj 120001 : 
iVan Waerbeke et al.l l2000; Kai ser et al.ll2000D . 

A direct measurement of the Hubble expansion rate 
at different redshifts can be obtained from the relative 
galaxy ages (RGA), i.e. the d ifferential age of galax - 
ies which are passively evolving ([Jimenez fcLoebl 120021 ) . 
This method was verified at low redshift with SDSS 
data ((Jimenez et al.l l2003f) , and has been used in the 
reconstruction of s calar field potential of dark energy 
(jSimon et alj l2005f ) and the constrain t on the cosmo- 
logical model firstlv (|Yi fc Zhand[2007f) . It is potentially 
competitive with other methods which probe cosmic ex- 
pansion history, such as SN la and BAO. Indeed, it mea- 
sures H(z) directly while SN la measures only the dis- 
tance and is related to H(z) by integration, hence the 
RGA meth od could be an even more sensitive probe in 
some cases dLin et al.ll2008l) . 

In this work, we focus on using WL data to constrain 
the equation of state of dark energy and the masses of 
neutrinos. With the current WL survey data, the SN 
la data, the baryon acoustic oscillation (BAO) data, 
and the Hubble parameter H(z) derived with the RGA 
method, we constrain the EOS w of the dark energy in a 
constant w model, the sum of the neutrino mass ^2 TO ^i 
and the other cosmological parameters. We also discuss 
the degeneracy of the parameters, especially for ^ m v 
and uu with the different data sets. The CMB data are 
not included in this work, as we wish to compare the re- 
sults from the weak lensing with that of the CMB data 
(WMAP 5-year measurement). We shall assume the ge- 
ometry of the Universe is flat in this work. 

The outline of this paper is as follows. In Sec. 2 we 
present our methods and the observational data sets. 
The results are given in Sec. 3. Finally we draw con- 
clusions in Sec. 4. 

2. METHODS AND OBSERVATIONAL DATA 

In this section, we briefly introduce the theoretical 
predictions of the observables. We employ the Markov 
Chain Monte Carlo (MCMC) technique to constrain the 
parameters. Eight MCMC chains are generated for each 
combination of our data sets, and after the convergence 
each chain contains about 100000 points to sample the 
probability distribution in the parameter space. Then, 
these chains are thinned and joined together, and there 
are about 1 2000 points left to be utilized to perform the 
constraints (|Gong fc Chedl2008h . 

2.1. Weak lensing 

As we know, the dark energy could affect the expansion 
rate and the large scale structure (LSS) of the Universe, 
from observation of these two aspects of cosmological 
evolution we can measure determine the property of dark 
energy very well. The WL provi des us a powerful probe 
for both aspects (jHuteren I2002T ). The massive neutri- 
nos could suppress the matter power spectrum on small 
scales, due to their free streaming, thus reducing the con- 



vergence power spectrum of the weak lensing, which is 
sensitive to the small scale matter distribution. Weak 
lensing is therefore a powerful measurement for both the 
dark energy and the massive neutrinos. 

For details of the WL tec h nique the readers can refer 
to iBartelmann fc Schneide"i1 (|2001h . We start with the 
power spectrum of the convergence k, which describes 
the strength of the lensing effect. Under Limber's ap- 
proximation, the convergence power spectrum P K can be 
related to the matter power spectrum P(fc, z) as 
c f Zs W 2 (z) 

r ^(,)r(z) 
~ I 3 Jo H E(z) A (M)dZ ' (1) 
where A 2 (fc, z) = -^jPik, z) is the dimensionless matter 
power spectrum, k = l/r(z) with I the multipole and r(z) 
the comoving distance defined as f Q -^fy, Ho is the 
Hubble constant, and E(z) is the expansion rate of the 
Universe 

e{z) = ^ = v/n m (i + z )3 + (i-n m )(i + 2)3(1+-), 

(2) 

where VL m is the matter density parameter and w is the 
EOS of dark energy. We give the detailed formula for 
the calculation of P K in the Appen dix. The shear cor- 
relation functions can be defined as (jMunshi et alJl2008t 
iDaniel et al1l2008l) 

£+(0) = ^ [°° dllP K (l)J (W), (3) 
2tt Jo 

£-(0) = ^ [°° dllP K (l)J 4 (W), (4) 
27r Jo 

e(0)=Ue) + if dtf'^P - 12^ 2 dtf'^fe) 

{.W-W^,^-^™, (6) 
where Jo and J4 are the zeroth and forth order Bessel 
functions of the first kind respectively. The shear cor- 
relation functions can then be compared with the mea- 
surem ents directly. In this work we use the £e(#) data 
from (Be njamin et al.l [2007) , which contains two wide 
sky survey: the Canada- France-Hawaii Telescope Legacy 
Surv ey of Wide field s (CFHTLS-Wide, sky coverage 22 
deg 2 . fFu et "aTll2008D and the Re d-Sequence Cluster Si ir- 
vey (RCS, sky coverage 53 deg 2 . iHoekstra et al.ll2003 ). 

The WL observations play the key role in our con- 
straint on the neutrino mass. We modify the code of 
calculating matter power spectrum to account for the 
suppression due to neutrino mass, the details of our cal- 
culation can be found in the Appendix. 

2.2. Type la Supernovae 

The SN la are widely used as standard candles 
to measure the l umin osity distance (|Riess etmil998l 
iPerlmutter et al.1 |T999) . The redshift-dependent lumi- 
nosity distances cIl(z) of the SN la are determined by 
the expansion history and geometry of the Universe. In a 
spatially flat Friedmann-Robertson- Walker (FRW) Uni- 
verse, the luminosity distance with redshift z is given by 

c (i + z) r- dz' ' 
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Then the distance modulus of the SN la can be written 

as 

IHh(z) = 5]0gd L (z) + 25. (8) 

We use here the SN la data recently published 
by the Supernova Cosmology Project (SCP) team 
(jKowalski et al.ll2008l ). This data set contains 307 SN la, 
selected from several current widely used SN la data sets 
inclu ding the Hubble Space Telescope (HST iRiess et al 



2004), SuperNova Legacy Survey (SNLS) (|Astier et al 



2006f ) and the Equation of Sta te: SupErNovae trace Cos - 
mic Expansion (ESSENCE, lWood-Vasev"et~aT1 12051 . 
They were re-analysed by the SCP team with the same 
procedure to get a consistent and high-quality "Union" 
SN la data set. 

2.3. Baryonic Acoustic Oscillation 

The acoustic oscillations in the plasma of the early 
Universe is imprinted on the matter power spectrum. 
This signatures in the large-scale clustering of galaxies 
yield additional cosmological tests. Using a large spec- 
troscopic sample of 46748 luminous red galaxies cover- 
ing 3816 square degrees out to z = 0.47 from the Sloan 
Digital Sky Survey (SDSS), Eisenstein et al. (2005) suc- 
cessfully found the acoustic peak in the matter power 
spectrum. The position of the feature can be described 
by the model-independent A-parameter 
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f zi dz' ' 
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(9) 



with Z\ = 0.35 the redshift at which the measurement 
is taken The value of the A parameter is measured as 
A = 0.469 ± 0.017 (jEisenstein et al.ll2005D . 

If the neutrinos are taken into account, the SDSS con- 
strain t on BAO can be approximated as (jGoobar et alJ 
2006) 



A = 0.469 



(—) 
V0.98/ 



-0.35 



(1 + 0.94/^ + 0.017, (10) 



where n s is the primordial power spectrum of fluctua- 
tions (see the Appendix), f v — Q, v /Q, m is the fraction of 
neutrinos relative to matter density. 

2.4. H(z) from relative galaxy ages 

The Hubble parameter H (z) is related with the differ- 
ential age of the Universe by 

1 dz 



H(z) = 



1 



At 



(11) 



so through the determination of dz/dt it can be measured 
directly. The later can be determined by using the differ- 
ential ages of passively evolving gala xies observed in the 
Gemini Deep Deep Survey (GDDS). iSimon et all (l2005h 
obtained a set of 9 H(z) measurement in the redshift 
range ~ 1.8 using this method. Various cosmologi- 
cal models were tested using this data set in the last 
few years (lYi fc Zhanel 120071: ISamushia fc Ratral 1 2006 : 



Wei fc Zhand 120071: iQiang et all 120071: iLin et al.l hoOS 
Carvalho et all l200ct lFigueroa et alJ I2008D . The con 



straint derived using this method is compatible and com- 
parable to that of SN la (ILin et al.ll2008t ICarvalho et~aTl 
[20081: lFigueroa et al.ll2008h . 



3. CONSTRAINTS ON COSMOLOGICAL PARAMETERS 

In this section, we first study the case of pure wCDM 
model, which contains a dark energy component with 
constant EOS w and the cold dark matter. The sepa- 
rate constraints for each data set are produced to illus- 
trate their power of constraints on the parameters and 
the degeneracy directions. Next we include the massive 
neutrinos component, first in the ACDM model (hence- 
forth ANCDM model), and then in the wCDM cosmol- 
ogy (wNCDM model). We give upper limits on vn v 
from the marginalized probability distribution function 
(PDF) , and discuss the degeneracy between the J2 m v 
and w. 

3.1. EOS of dark energy in wCDM model 



* -1 - 




Fig. 1. — The contour maps of Cl m -w for the WL (red solid), SN 
la (blue dash-dotted), RCA (magenta long dashed), BAO (green 
short dashed) and the combination constraint (thick black solid). 
For each data set the lcr (68.3%) and 2cr (95.5%) confidence levels 
(C.L.) are shown respectively. 

The la and 2a contours of the f2 m -u> plane are shown 
in FigQ] We find that the WL data is sensitive to the 
matter density parameter fl m , but less effective in con- 
straining the dark energy EOS w. Also, the best-fit 
r2 m is nearly independent of the w for the WL data, 
this is consistent with what were found in other works 
([Hoekstra et alj|2006t iBeniamin et~aT]|2007l) . so the de- 
generacy direction of the WL is different from other cos- 
mological observations. It is also shown that the degen- 
eracy directions of the SN la data and the RGA data 
are almost the same, which again agr ees with the con - 
clusions reached in th e recent works of [Lin et al.1 (|2008h : 
ICarvalho etldl (|2008h : lFigueroa et"aH (|2008h . The con- 
straint from the SN la data is much tighter than that 
of the RGA data due to the much larger sample of the 
former. If we combine all of the data sets, the results 
are significantly improved. We measure from the global 



fit of WL+SN Ia+RGA+BAO that, VL m = 0.28 



+0.04 
-0.03 



and 



w = — l.OOlto i2) which are consistent with the recent re- 
ported results from WMAP 5-year data (iKomatsu et al.l 
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Fig. 2.— The marginalized PDF of J2 m v in ANCDM model. 
The 2cr(95.5%) C.L. is also shown by the blue dashed line. 



2008). This result show that the cosmological constant 
is an excellent candidate of dark energy. 

3.2. Constraints on neutrino mass 

Now we investigate the constraints on neutrino 
mass for the employed data sets. We use the WL, 
the WL+RGA+BAO, the WL+SN Ia+BAO and the 
WL+SN Ia+RGA+BAO data set combinations to per- 
form the constraints, so that we could compare the power 
of parameter constraint for the different data set combi- 
nations. 

Due to parameter degeneracy, with WL only there 
is very little constraining power. With the addition of 
BAO and RGA data, a 2a limit of 0.6 eV can be ob- 
tained. With SN la instead of RGA the constraint is 
better (0.4eV). If all (WL, BAO, SN la and RGA) data 
are combined to make the constraint, it is further slightly 
improved. It is also notable that the peak of the PDF 
distribution is not at eV but at about O.leV for the 
last three data set combinations. 

In Fig. [3l we show the contour map of the £l m -w plane 
in the wNCDM model. As can be seen, the results do not 
change much compared with that of the wCDM model. 
The constraint of w for WL+RGA+BAO is looser than 
that for WL+SN Ia+BAO, but these two combined data 
sets yield similar constraint on Q m . This comparison ac- 
tually reflects the constraining ability of the SN la data 
and the RGA data. As we shall also see in other results 
presented below, except for w, the WL+RGA+BAO al- 
most have about the same constrain ing power as the 
WL+SN Ia+BAO data (|Lin et al.ll2008ft . even though the 
number data points in the RGA data set is much less 
than that of the SN la data set. 

The contour maps for w„ are shown in Fig. [3J 
The constraint is not good if only the WL data are em- 
ployed. After including the other data sets, the con- 
straints of ^2 m„ are improved much, because the SN la, 
RGA and BAO data could remarkably improve the con- 
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Fig. 3. — The contour maps of Q. m -w for the different data sets. 
The 1o-(68.3%),2(t(95.5%) and 3<x(99.7%) C.L. are marked by red 
solid, green dashed and blue dotted lines respectively. 
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Fig. 4. — The contour maps of m-^m,, for the different data 
sets. The 1ct(68.3%), 2cr(95.5%) and 3<r(99.7%) C.L. are marked 
by red solid, green dashed and blue dotted lines respectively. 

straints on the other cosmological parameters such as Q m 
and w. We can also see that Yl m u tends to be greater 
when the w becomes more negative; however, there is 
no strong degeneracy between w and ^2 u t . The reason 
of this may be, on one hand, the current WL data are 
not yet accurate enough to indicate such relations; on 
the other hand, the P K , which is used to constrain the 
parameters for the WL data, is the integral of the mat- 
ter power spectrum P{k,z), so it is not as powerful as 
P(k, z) t o reflect the influence of w on the fo rmation of 
the LSS (|Hannestadll2005t ISpergel et al.ll2007ft . 

Finally, we marginalize over the other parameters and 
get the PDF of m v shown in FigJSJ the 2a C.L. is also 
shown as an vertical line. We find 'Y^m v < 0.8 eV at 
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Fig. 5. — The marginalized PDF of the ^2m v for the different 
data sets. The 2<r(95.5%) C.L. is also shown by the blue dashed 
line. 

95.5% C.L. for the WL+SN Ia+BAO and the WL+SN 
Ia+RGA+BAO data. This result is comparable with 
the recently reported results TO f — 0.66 eV from a 
combined analysis of WMAP 5-year data to gether with 
SN la and BAO data (jKomatsu et alJl2008f ). Moreover, 
as can be seen from this figure the constraint range of 
^2 m v is fairly large if we just use the current WL data, 
but the PDF for the WL+RGA+BAO data set is almost 
identical to that for the WL+SN Ia+BAO, as we have 
discussed earlier. 



WL along does not provide a constraint, but if one com- 
bine WL with SN, BAO or RGA data, the neutrino 
mass is well constraint. The 2a limit with all of these 
observations together is about 0.4 eV, and the peak of 
the fit is around 0.1 eV. We also considered the wCDM 
model with massive neutrinos. However, we find that 
w = 1 (cosmological constant) is still favored in this 
model, and a weak degeneracy is found between w and 
Y^ m v- It shows that WL can indeed break the degen- 
eracy between w and m u effectively (|Hannestad et all 
2006). After marginalizing over the other parameters 
we obtain the PDF for ^m„. The upper limit of the 
neutrino mass is obtained as Yi m v — 0-8 eV at 2cr, 
which is comparable with t he recent combined a nalysis of 
CMB+SN Ia+BAO data (jKomatsu et al.ll2008l ). More- 
over, except for w, the constraint ability for the WL+SN 
Ia+BAO and the WL+RGA+BAO data set is nearly 
the same, which means that RGA data can play a simi- 
lar role in cosmological study as SN la (|Lin et al.ll2008t 
ICarvalho et~alll2008t iFigueroa et al.ll200iT " the p oint of 
which is clearly demonstrated bv lLin et alj <|2008f> . 

Although the capability of the current WL data 
on cosmological parameter constraints is still not very 
good, and our constraints on w and m v ar e not as 
strong as the other works (ISeliak et al.l l2005alfbl. 120061: 
Kristiansen et al.1 l2007: Ko matsu et al.ll2008t iFogli et al.l 
20081 ). it shows that WL is a good and powerful comple- 
mentary for the other observations, and plays a more and 
more important role for the cosmological study. When 
the next generation WL surveys, such as the Super- 
Nova/ Acceleration Probe 1 (SNAP) and the Large Syn- 
optic Survey Telescope 2 (LSST), are put to work, the 
WL observation would be more accurate and become an 
indispensable measurement for the cosmology study. 



4. CONCLUSION AND DISCUSSION 

In this work we use the WL survey data, combined 
with SN la, RGA and BAO data to constrain the EOS of 
dark energy w and the sum of the neutrinos mass ^ m„. 
The MCMC method is employed to give the PDF of the 
parameters. 

We first study the wCDM model without the neutri- 
nos, and obtain the constraint on fl m -w for different cos- 
mological observations. We find that the constraint for 
w is weak if only the WL data are used, but after in- 
clusion of the other observational data, the constraints 
are improved remarkably, and the cosmological constant 
(w = — 1) is favored. 

Then we considered constraints on the neutrino 
masses. In the ACDM model with massive neutrinos, 



1 http://snap.lbl.gov 

2 http://www.lsst.org 
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APPENDIX 

CALCULATION OF LENSING POWER SPECTRUM 
The weight function W(z) in Eq.(TJJ is written as (H utererll2002f ) 



2 V c 

with function 



W(z) = U^) n m f(z)(l + z), (Al) 



/(--) = /(:)/ ^'l^ n(z')dz\ (A2) 



in which n(z) is th e normalized number d ensity distribution (i.e., J n(z)dz — 1) of source galaxies. The n(z) we use 
here is proposed in lBenjamin et ail (|2007l ). which takes the form as 



n(z) = N 



z b + c 



where a, b and c are free parameters, and N is a normalizing factor, 



dz'- 



The linear matter power spectrum A^(fc, z) is parameterized as 



(A3) 



(A4) 



A£(k, z) = Ak n ° +6 T\k)D\z), 



(A5) 

where D(z) = g(z)/(l + z)g(0) is the linear growth factor, T(k) is the transfer function, A is the normalization factor 
and n s is the primordial fluctuation spectrum. Hereafter we use a Harrison-Zel'dovich spec trum n s = 1. For the 
ACDM model (w — —1), the relative growth factor g(z) is found to be well approximately by (j Carroll et al.lll992l ) 



with 



9a(z) 



(5/2)n m (z) 



nU\z) - n A ( z ) + (i + n m {z)/2)(i + n A (*)/70) 



fim(l + zf 



SI, 



n m (i + z)3 + o A J1W n m (i + zf + n A 

For the transfer function, we adopt the fitting result of iBardeen et al.l ([1986) for an adiabatic ACDM model 



T A (q) = 



ln(l + 2.34g) 
2.Mq 



[1 + 3.89q + (16. Iqf + (5A6q) 3 + (6.71g) 4 ]" 1/4 , 



(A6) 



(A7) 



(A8) 



where q = k/hT, and h = Hq/(100 kms^ 1 Mpc _1 ),and T — fi m /iexp[— Qf,(l + V2h/Q m )] is the shape parameter 
with baryon density Ob . fib = 0.0462 is adopted according to the recent analysis of CMB, SN la and BAO data 
(|Komatsu et aLlfeOOSf) . For the extension of the growth fact or and transfer f unction from ACDM model to any dark 
energy case with constant EOS w, we use the fitting form of iMa et"aT1 (|1999f ). 



9Q 
t 



9x{-w)\ 
-(0.255- 
-(0.366- 

Q m )(l + z) 3w ] 



0.305w + 0.0027/w)[l - fi™ (z)] 
0.266w-0.07/w) In n™(z), 

This fitting formula is accurate to 2% for 0.2 < Vt r 



where Cl™(z) = 
-1 < w < -0.2. 

For the non-linear power spectrum, we adopt the formula given bv lPeacock fc Doddsl (|1996h . 



(A9) 

(A10) 
< 1 and 



Aj, L (k NL ) = f NL [A 2 L (k L )} 
k L = [l 

1 + B(3x 



A^ikNL)}' 1 / 3 ^^ 



/nl(x)=x 



(Ax) 



a0 



1 + [(Ax)<*g 3 (z)/(VxV 2 )]P 



VP 



(All) 



The parameters in the non-linear function f^L are 

A = 0.428(l 
5 = 0.226(1 
a = 3.310(l 
^ = 0.862(1 
V = 11.55fl- 



n s /3)- a947 
n./S)- 1 - 778 
n s /3)-°- 244 
n s /3)- a287 
n s /3)-°- 423 



which are fitted to numerical simulation results. 

Actually, the non-linear matter power spe ctrum we emplo y can be seen as the Peacock-Dodds (PD96) fitting formula 
with a modified growth factor propose d bv (IMa et al.lll999f). There is not yet a reliable analytical fit of the non-linear 
power spectrum for the wCDM model (|Hoekstra et al.ll2006f ). so now it is usually obtained from fitting to the N-body 
simulations, the applicable range in the parameter space i s always small. F o r inst ance, an accurate fitting non-linear 
matter power spectrum with varying w was proposed by McDon ald et al.l (|2006l ). but unfortunately, its applicable 
range was Q, m G [0. 211, 0.351] . However, We have tested our power spectrum and find it matches well with the 
simulation results of iMcDonald et all (|2006D , especially for w < —0.5 (see Figj6]). 



0.1 



1 

k [hMpc- 1 ] 
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Fig. 6. — The ratio of the mCDM matter power spectrum P m to the ACDM matter power spectrum Pa as a function of k. The P w 
is calculated by our non-linear power spectrum code. The dash-dotted, dashed and dotted lines represent w = —0.5,— 0.75 and —1. 5 
respectively, and we set £l m = 0.281 which i s to match the paramet er value for the black solid line in Fig. 1 of McDonald ct al. (2006). 
We find our results match well with that of McDonald ct al. (2006) when w < —0.5. Even for w = —0.5, our P w /Pa is about 1.6 at 
k = 10hMpc~ 1 while it is about 1.35 in[McDonald et al. (2006), that still does not deviate much. 

If there is a fraction of massive neutrinos in the matter components of the Universe, the growth of the structure is 
suppressed by the free streaming of neutrinos. The transition scale is the horizon scale when the neutrinos become 

non-relativistic: k nr » 0.026 (^) 1/2 fi^/iMpc -1 , below which (|Hu fc Eisensteinlll998fl 

AP L n Q v 



Pl fi, 



(A12) 



where £l v — J2 m ly /(93.2eV h 2 ) is the neutrino matter density (jGoobar et al.ll2006f) . This approximation is shown to 
be accurate for the lin ear theory regime with relative large scales k < 0.2h Mpc -1 and small neutrino fraction f v 
(Brandbv ge et al] 12008). Hence, When we consider the neutrino component the linear matter power spectrum should 
be modified as P£ — Pl + AP^, (k > k nr ). We add this modification into our non-linear matter power spectrum code 
to constrain the sum of the neutrino mass ^m v . 
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